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Abstract: The analysis of the spatiotemporal trends of temperature and precipitation are important for 
the policy guidelines of future development and sustainable management of a country. This study 
presents spatiotemporal variations in temperature and precipitation indices to elucidate climatic change 
over Iran during last five decades i.e., 1961 -2010 by examining data obtained from30 synoptic stations 
across the country. Results indicate significant and spatially coherent rising trends in temperature 
indices. From the results it is clear that the majority of the trends in the annual means of daily 
maximum temperature, daily minimum temperature, tropical nights, warm nights and warm spell 
duration index showed strongly positive tendency during the last decades, while the increasing trends in 
the Central region of Iran was stronger than this in the other areas. Decreasing trends have been found 
in the diurnal temperature range, Cold spell duration index and Cool nights were found over most 
regions of Iran and positive trends were found for Growing season length and Warm days over other 
regions of the country. Cooling trends were found for the southern half of the Zagros Mountain and 
Shahrekord station. Among the precipitation indices we found positive trends in daily intensity index, 
heavy precipitation and extreme wet days over most of the country and a significant increase in the 
heavy precipitation series was observed mostly in Southwestern regions of Iran, as well as along the 
coasts of the Persian Gulf. We can conclude that the average amount of wet day precipitation has 
significantly increased across Iran. Overall, these trends were statistically significant mostly over the 
central areas of the country and there was strong climatic variability in the rest of country, especially 
arid and semi-arid regions. The northwest region was only region with weak anomalies in studies 

indices. 
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1. Introduction 

The subject of climate change is gaining 
importance due to increasing frequency, duration and 
intensity of extreme climate events, such as drought, 
flood, heat wave and cold waves. Extreme heat events 
includes record breaking summer heat waves of 2003 
and 2010 hampered various aspects of life due to 
intensity, duration and frequency of high temperature 
in different parts of the world (Garcfa-Herrera et al., 
2010; Barriopedro et al.,2011; 

Coumou and Rahmstorf, 2012). Numerous studies 
using a long term dataset showed that global warming, 
as a result of greenhouse effect, has very likely 
increased temperature variability in Europe, as a 
result more frequent and severe warm extremes were 
recorded while there was an attenuation of cold 
events ( Keggenhoff et al., 2006; Fischer and Schar, 
2010; Trenberth and Fasullo, 2012; Krueger et al., 
2015). The study of these events and their frequency 
is very important in order to understand causes and to 
decipher certain signs or particular pattern of their 
occurrence. Moreover, this kind of research is 



important for the assessment of their impacts on 
human beings, agriculture and industries. 

IPCC (2007) described 0.74 °C rise in global 
mean surface temperature froml906 to 2005. This 
increase in temperature was more significant 
duringl910 to 1945 and 1970 to 2005. Other studies 
found consistent increase in temperature since 1980s. 
The analysis of temperature during last 50 years 
showed that the warming occurred at the rate of 
0.13°C per decade, which is almost double than the 
centennial average of warming at the rate of 0.07°C 
per decade over the last 100 years. Moreover it was 
observed that, between 1950s and 1990s northern 
hemisphere was warmer than any other 50-year 
period in the last five centuries. It is reported that a 
reduction in consecutive dry days (CDD) and warm 
nights were observed. Asymmetric rise in temperature 
is the prominent feature, where minimum temperature 
is increasing faster than the maximum temperature, 
resulting in the decreased diurnal temperature range 
(DTR) at the rate of -0.07°C per decade since 1950s, 
(Boccolari and Malmusi, 2012). 
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Most of the studies had analyzed the long term 
climatology data with the focus on mean values, 
neglecting extreme climatic events. This was mainly 
due to unavailability of high quality meteorological 
data which is usually prerequisite for monitoring, 
detection and attribution of anomalies in climate 
extreme (Zhang et al., 2005). Recently Rahimzade et 
al., (2009) analyzed meteorological data of Iran using 
27 different indices of temperature and precipitation, 
as indicator of climatic anomalies to determine recent 
climate change. They found fairly negative trends for 
low temperature or winter indices like cool days and 
nights, frost and ice days, and DTR. However 
positive trends were found for warm days, summer 
days, and tropical nights over most region of Iran. 

In addition, the analysis conducted by different 
researchers in different countries may not perfectly 
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merge to form a global map as they might have 
investigated different indices (Bonsai et al., 2007) and 
potentially using different range of thresholds. 
However, anomalies and trends in climatic extremes 
and extreme events have shown more sensitivity to 
climate change as compared to the means values. 
Extreme events have more serious implications than 
average trends therefore they have received more 
attention from scientific community (You et al, 2010a, 
2010b). Throughout the 20' century anomalies in 
climatic extremes had been significant across the 
globe. Among these extremes rising temperature was 
pronounced and precipitation indices showed a 
tendency toward wetter conditions particularly in the 
last few decades. (Frich et al., 2002; Alexander et al., 
2006; Coesar et al., 2006). 
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Fig. 1. Spatial distribution of the synoptic stations in Islamic Republic of Iran map. 
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Fig. 2. Five regions including: Azarbaijan (A), Zagros (B) and Alborz (C) mountain ranges, Caspian Sea cost 

(D) and central part (E) of Islamic Republic of Iran. 



Both temperature and precipitation extremes have 
teleconnections with each other's (Krishnamurti et al., 
2015). Pakistan flood and Russian heatwave of 2010 
are recent example of teleconnection between 
hydrometeorological extremes (Lau and Kim, 2012) 
Studies have shown a significant widespread changes 
in temperature extremes are associated with global 
warming. Changes in precipitation extremes exhibit 
much less spatial coherence than those of temperature 
extremes, making it different to detect regional trends 
(Haylock et al, 2006; Klein Tank et al, 2006; 
Trenberth and Fasullo, 2012). 

This study was aimed to provide a better 
understanding of the anomalies and changing trends 
in duration, frequency and intensity of climatic 
extremes by examining temperature and precipitation 
indices over Iran during 1961-2010. 

2. Material and Methods 
2.1. Data 

The study area i.e., Iran, lies approximately 
between 25 °N and 40 °N in latitude and between 
44 °E and 64 °E in longitude (Fig.l). Based on the 



Koppen climate classification, most parts of Iran are 
categorized under arid (BW) and semi-arid (BS) 
climates. Alborz and Zagros are the important 
mountains of Iran, which play an important role in 
non- uniform spatiotemporal distribution of 
temperature and precipitation in Iran (Dinpashoh, 
2006). 

Examination of climate changes needs long and 
high quality records of climatic variables. In the 
present study, dataset of daily maximum, minimum 
and mean air temperatures and precipitation (P) for 
the period 1961-2010 from 30 synoptic stations from 
different geographic locations of Iran (Fig.l) were 
collected from the Islamic Republic of Iran 
Meteorological Organization (IRIMO) and were 
analyzed. Homogeneity of the dataset was assessed 
and approved by IRIMO previously. Most of the 
regions were covered by the corresponding data and 
the geographical location of the stations. 

2.2. Regional categorization 

To provide an overall picture of climate variation 
in the Iran, we also computed regional averages for 
every index. We divided total area of Iran into five 
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climatic regions based on the climatic variation and 
geographic locations i.e., distance from sea, 
mountains and central land. These regions includes 
mountain ranges Azerbaijan (A), Zagros (B) and 
Alborz (C), coastal areas of Caspian Sea (D) and 
central part of country (E) (Fig. 2). The regional 
averaged indices are computed as the mean of the 
indices at individual stations relative to their 
climatology. 



Tablel. List of 29 ETCCDMI Core Indices for the analysis of extreme precipitation and temperature. 



ID 


Indicator Name 


Definitions 


Unit 


TXMean 


Annual maximum temperature 


Annual mean of daily maximum temperature 


°C 


TNMean 


Annual minimum temperature 


Annual mean of daily minimum temperature 


°c 


FDO 


frost days 


Frost days: Annual count when TN (daily minimum) < 0 C 


Day 


SU25 


summer days 


Summer days: Annual count when TX (daily maximum) >25 C 


Day 


IDO 


jpp Havs 


Tpp Havs' Annual rnunt whpn T 1 "^ fHailv maximum^ ^0 C 


Day 


TR20 


tropical nights 


Tropical nights: Annual count when TN (daily minimum) >20 C 


Day 


GSL 


growing season length 


Growing season length: Annual (1 Jan-31 Dec in Northern Hemisphere) 


Day 






count between fir st span of at least 6 days with daily mean temperature 








T>5°C. 


C 


TXx 


max Tmax 


Max Tmax: Monthly maximum value of daily maximum temperature 


TXn 


min Tmax 


Min Tmax: Monthly minimum value of daily maximum temperature 


c 


TNx 


max Tmin 


Max Tmin: Monthly maximum value of daily minimum temperature 


c 


TNn 


min Tmin 


Min Tmin: Monthly minimum value of daily minimum temperature 


c 


TNlOp 


cool nights 


Cool nights: Percentage of days when TN < 10th percentile 


Days 


TXlOp 


cool days 


Cool days: Percentage of days when TX < 10th percentile 


Days 


TN90p 


warm nights 


Warm nights: Percentage of days when TN > 90th percentile 


Days 


TX90p 


warm days 


Warm days: Percentage of days when TX > 90th percentile 


Days 


WMJ1 


warm spell duration indicator 


Warm spell duration indicator: Annual count of days with at least 6 


Days 






consecutive days when TX > 90th percentile 




CSDI 


cold spell duration indicator 


/""I U 11 J t-' J ' d. A 1 i P J "i-l i-l i- S~ 

Cold spell duration indicator: Annual count of days with at least 6 


Days 






consecutive days when IN < 10th percentile 


0 

c 


DTR 


diurnal temperature range 


Diurnal temperature range: Monthly mean difference between TX and TN 


CDD 


consecutive dry days 


Consecutive dry days: Maximum number of consecutive days with 


Days 






precipitation < 1 mm 




CWD 


consecutive wet days 


Consecutive wet days: Maximum number of consecutive days with 


Days 






precipitation >lmm 




PRCPTOT 


annual total wet day 


Wet days precipitation: Annual total precipitation in wet days with 


Mm 




precipitation 


precipitation >lmm 




RlOmm 


number of heavy 


Heavy precipitation days: Number of days with precipitation >10 mm 


Days 




precipitation days 






R20mm 


number of heavy precipitation 


Very heavy precipitation days: number of days with precipitation >20 mm 


Days 




days 






R25mm 


number of heavy precipitation 


Very heavy precipitation days: Number of days with precipitation >25 mm 


Days 




days 






RxlDAY 


number of heavy precipitation 


Max 1-day precipitation: Monthly maximum 1 day precipitation 


Mm 




days 






Rx5DAY 


max 5 days precipitation 


Max 5-day precipitation: Monthly maximum consecutive 5 days precipitation 


Mm 




amount 






SDII 


simple daily intensity index 


Simple daily intensity index: Annual total precipitation on wet days divided 


mm/ 






by number of these days with precipitation >lmm 


day 


R95p 


very wet days 


Very wet days: Fraction of annual total precipitation exceeding the 1961- 


% 






1990; 95th percentile 




R99p 


extremely wet days 


Extremely wet days: Fraction of annual total precipitation exceeding the 


% 






1961-1990; 99th percentile 
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2.3. Data analysis techniques 

In the study a linear trend was computed from the 
indices series using a least squares fitting. The Mann- 
Kendall statistic tests were applied to detect and 
describe significant trends in the climatic variables. In 
addition, the magnitude of the decadal trends was 
derived from the Durbin-Watson test. 
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The brief descriptions of the used statistical methods 
are as follows: 

2.3.1. Mann- Kendall test 

The Mann- Kendall test is one of the widely used 
non- parametric tests to detect significant trends in 
hydrological and meteorological time series data 
(Modarres and da Silva, 2007; Partal and Kahya, 
2006). The test compares the relative magnitudes of 
sample data rather than the data values themselves 
(Gilbert, 1987). Mann (1945) originally used this test 
and Kendall (1975) subsequently derived the test 
statistic distribution. The test is suitable for data that 
do not follow a normal distribution, and supports 
multiple observations per time series (Kampata et al., 
2008). The Mann- Kendall test is given as: 



■5=^^ sign {%j - 



fc=U=te-rl 



(2) 



(1) 



+1 if fa -xj >0 
-1 if (x } - x k ) < 0 



Var© = 

(3) 



IS 



Where n is the number of data points, t is the 
number of ties for i value and m is the number of tied 
values. Then equations (1) and (3) were used to 
compute the test statistic Z from the following 
equation: 

S-t _ 
_ ifS > 0 

: , ={o ifs=m 

.JVARift (4) 

A positive value of Z indicates an increasing trend 
and a negative value indicates a decreasing trend. The 
null hypothesis Ho that there were no trends in the 
records is either accepted or rejected depending if the 
computed Z statistics is less than or more than the 
critical value of Z-statistics obtained from the normal 
distribution table at the 5% significance level 
(Kampata et al., 2008). 



2.3.2. Durbin- Watson test 

This test is based on the first order autoregressive 
error model e t = (pe,_i+w, Where cp is the 
autocorrelation parameter and the u t , are 
independently normally distributed with zero mean 
and variance a 2 . When one is concerned with positive 
autocorrelation the alternatives are given as follows: 

(5) 

Here Z/oimplies that error terms are uncorrected 
or negatively correlated, while Hi implies that they are 
positively auto-correlated. This test is based on the 
difference between adjacent residuals e,-e,_iand is 
given by 



d = 



(6) 

Where e, is the regression residual for period t, and 
n is the number of time periods used in fitting the 
regression model (Kanji, 2006). 

2.3.3. Selected indices 

The ETCCDMI recommended a total of 29 core 
indices with primary focuses on extremes to be 
derived from station daily data (Table 1). We 
compute these indices using RClimDex, an R-based 
software package developed at the Climate Research 
Branch of Meteorological Service of Canada on 
behalf of the ETCCDMI. 

3. Results and Discussion 

3.1. Extreme Temperature 

The results highlighted that the annual means of 
daily maximum temperature (TX), minimum 
temperature (TN), summer days (SU25), tropical 
nights (TR20), highest daily minimum temperature 
(TNx), highest daily maximum temperature (TXx), 
lowest daily maximum temperature (TXn), lowest 
daily minimum temperature (TNn), warm days, warm 
nights, duration of warm spell and growing season 
length showed increasing trend in most parts of the 
country (Fig.4 and Fig. 6). 

The positive trend of annual means of TX and TN 
especially since 1970s is an important indicator of 
warming trend in the country (Fig. 4). These findings 
are in good agreement with the results obtained from 
other studies in the Iran (Tabari et al., 201 1). 

Significant positive trends of TR20 and summer 
days (SU25) for mountainous and arid and semi-arid 
climatic regions were recognized. The largest positive 
trends of summer days (SU25) about 1.6 and 1.5 day 
per decade occurred in Babolsar (Fig.3A) and 
Sabzevar stations. But the largest negative trend for 
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this index of -2.09 and -1.01 days per decade decade) (Fig.3B) and Ghazvin (-1.3 days per decade) 



occurred in Anzali and Saghez stations. 

These findings support results of Zhang et al. 
(2005). 21 stations in the southern and central regions 
have experienced positive trends for tropical nights 
(TR20), especially after 1980s. Ahvaz (+1.5 days per 



showed higher anomalies. 

The rate of the positive trends in TXx (Fig. 4) and 
TNx (Fig.6) were obtained in this study was weaker 
than that reported by Zhang and et al (2005). Most of 
the studied stations i.e., 18 of 30 stations, 18 of 30 
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Fig. 3. Time series of temperature extreme of five different stations as representatives of different geographic 
locations. 



stations, showed a positive trends of TNx. The largest 
positive trend of 1.6 C per decade was observed in 
Anzali followed by Shahrood, The statistical tests 
detected a warming trend in the TXx at more stations 
such as in southern region. But negative trends were 
overall occurred in Zagros Mountain. Generally Iran 
experienced a positive trend for TXn (Fig.4) and TNn 
(Fig.6). Especially, most of the stations show positive 
trends for TNn, which lies approximately between 29°N 
and 33 °N and over 36 °N in latitude. The strongest 
positive trends for these daily temperature extremes 
have been found in Ahvaz and Babolsar with trends of 



2.1 and 1.9 C per decade, respectively. The significant 
negative slopes of linear trends for this index were -1.5 
C per decade in Bandar Abbas and -2.3 C per decade in 
Birjand. In regions such as the Zagros Mountain, the 
negative trends obtained are not so strong. 

Any substantial differences were not found between 
the TN90p and TX90p indices and many stations show 
significant positive trends (fig.5). These results are more 
consistent with the findings of Zhang et al. (2005) and 
there is a strong evidence for globally increasing 
tendency of warm nights both in terms of their number 
and intensity (Christidis et al., 2005). 



6 



Parak et al., 2015. 2:11. Journal of Environmental and Agricultural Sciences (ISSN: 2313-8629) 



(A) 



_ _ 



• 



TMAXmean 



- »<•• »n 
«•» *m 



— • 



• 



..." " 



(B) 



V 



♦ 



TMINmeao 



£3 «rr 



(C) 



FDO 



MM 



_ » 



(E) 




IDO 



«« OHO * 




(D) 



SOS 



^ *M 

■ — 

* .. « . . ... 





— • . 



(F) 




(O 



TXn 

Q • >■ • - 

Q •«-•» 

• « Ml 



(H) 



TXx 



_ • 



• M~ 



Figure 4. Trends in averages of TXmean, TNmean, FDO, SU25, IDO, TR20, TXn, TXx over Iran. 
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Fig.5. Comparison between time series of GSL for Zanjan (A), Shiraz (B) and CSDI for Shahrekord (C), 
Saghez (D). 

The patterns for both daily temperature extremes 
i.e., TN and TX at 90 th percentile (TN90p and TX90p) 
indices are quite similar and positive trends at most of 
the weather stations are statistically significant. 
Trends in the TN90p and TX90p were larger than 
those in the TNI Op and TXlOp. This Similar trends 
of these indices were also reported for the Middle 
East regions (Zhang et al., 2005). A similar pattern 
for TR20 was also found over the country (fig.3). 



The highest numbers of TX90p and TN90p have 
also occurred after 1980s, coinciding very well with 
the result of Zhang et al. (2005). We cannot found 
significant negative trends for TN90p over the whole 
country. More than 22 stations have experienced 
positive trends, the strongest being in Bam with 0.99 
percent in decade (Fig.3C). The result shows the 
decreasing and increasing trends of warm days for the 
Esfahan and Shahrekord stations, respectively. The 
number of warm days at the Shahrekord station had a 
decrease over the study period of approximately -1.33 
percent per decade while the Esfahan station had an 
increase of 1.01 percent per decade. 

The positive trends for the warm spell duration 
index (WSDI) were found in most parts of the 
country however this was prominent for Yazd with 



1.37 days per decade (Fig.3D). The rate of the 
positive trend in WSDI obtained for this study is very 
smaller than that (9 day per decade) reported by 
Rahimzade et al (2009). Negative but weak and non- 
significant trends of WSDI were observed for Anzali 
and Shahrekord stations. The trends of growing 
season length (GSL) were not meaningful for the 
weather stations located at the southern parts of the 
country. However negative trend of GSL was 
significant for Zanjan (Fig. 5A) and Saghez. 

On the other hand, Frost days, Ice days, Cold 
nights, Cold days, Cold spell duration and Diurnal 
temperature range decreasing in most parts of the 
country. Excluding of the southern region, negative 
trends in FD according to global decrease trend 
(Zhang et al. 2005; Frich et al. 2002) were observed 
from 1980s for other areas of the country. In some 
regions such as Southern Zagros (west) and Central 
and Eastern Alborz (north), negative trends were 
significant for the number of FD. Large significant 
negative trends in this indices were found in Yazd, 
Kerman, Sabzevar and Shahroud stations located in 
the arid regions of the country imply that frost days 
have decreased by about respectively 1.5, 1.5, 1.2 and 
1.1 days decade" 1 
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The rate of the negative trend in FD obtained in 
this study for Sabzevar station (-1.2 day per decade) 
is very smaller than that (- 9.9 day per decade) 
reported by Rahimzade et al. (2009). On the contrary, 
the magnitude of the positive trend in FDO found in 
this study only for Shahrekord station is 1.2 day per 
decade. The mountainous region of the country more 
experienced a maximum temperature below 0°C. In 
the mountainous region in the west and north i.e., 
Zagros and Alborz, contrasting trends for the number 
of FDO were observed, positive for Zargos, while 
negative for Alborz (Fig.4). 

Countrywide significant negative trends for 
TNI Op and TXlOp indices were observed. This result 
is more consistent with the studies by Zhang et al. 
(2005) and Rahimzade et al. (2009). Negative trend 
rates in cooling night lay in the range of -1.98 percent 
decade" 1 (Ahvaz station) to -1.4 percent decade" 1 
(Bushehr station), while decreasing trend in cooling 
day was in the range of 1.28 percent decade" 1 at 



Kermanshah station. On the contrary, the magnitude 
of the positive trend found in this study for 
Shahrekord station (0.93 percent decade A ) and 
Saghez station (0.8 percent decade " ). 

In case of cold spell duration indicator (CSDI) 
most of the stations showed negative trends for and 
more of them are statistically significant (Fig. 6). 
Saghez showed a negative trend for WSDI and 
significant positive trend for CSDI about 1.4 days per 
decade (Fig.5D). Although the slope of CSDI trend at 
Sharekord is about 1.13 days per decade (Fig.5C), 
while this station showed a negative trend for Warm 
spell duration indicator (WSDI). Most of the stations 
show negative trends for diurnal temperature range 
(DTR) due to decrease for monthly mean difference 
between daily temperature extremes i.e., TX and TN. 
The largest negative trend of -1.1, -1.03 and -0.95°C 
per decade occurred in Ramsar, Shiraz and Bam 
(Fig.7), however some of the stations show non- 
significant positive trends for this index. 
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Figure 7. Time series of negative trends for four temperature extremes in different stations. 



The increasing trends in temperature extremes i.e., 
Tmax and Tmin were found in this study are in 
accordance with the results of Tabari et al. (2011). 
But trends were overall stronger in the minimum 
temperature. Also, Marked significant trends for 
indices like Warm spell duration (WSDI), Cold spell 
duration (CSDI) and Growing season length (GSL) 
were found over most regions of Iran are in 
accordance with the results of Zhang et al. (2005) 
3.2. Extreme precipitations 

Positive trends were found for wet-days 
precipitation (PRCPTOT), extremely wet days 
(R99p), very wet days (R95p) and simple daily 
intensity index (SDII), and over most regions of the 
country. Among the studied 30 stations majority of 
stations showed positive trends and the 10 stations 
showed significant trends for SDII. The results 
highlighted that annually occurrence of very wet days 
(R95p) during the studied period in stations located at 
the Caspian Sea shore. 

Positive trends in R95p at stations from arid and 
semi-arid climate regions of Iran, in addition verify 
the presence of the more extreme precipitation events. 



The largest significant positive trends occurred in 
Esfahan (1.7 percent decade 1 ) and Tehran (2.57 
percent decade" 1 ) respectively (Fig.9 A, B). The 
western regions of the country experienced a positive 
trend for extremely wet days (R99p), but significant 
trend only occurred in Abadan. Also, negative trends 
in R99P were found at the majority of the stations in 
the western, eastern and arid and semi-arid regions. 
The strongest negative trends were also found in Bam 
and Tabriz. The amounts of annual total precipitation 
and annual total precipitation from wet days 
(PRCPTOT) are very similar and substantial change 
of such wet days occurred in 1970. 

The significant negative trends for this index were 
found in the northwest and southeast regions of the 
country and southeast of the Caspian Sea, and 
significant positive trends in the north and southwest 
regions. In addition, the strongest positive trends in 
Babolsar and Abadan with trends of 2.2 and 2.1 mm 
per decade. On the contrary, the largest negative trend 
found in this study for Tabriz station (1.85 mm per 
decade). 
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Fig. 8. Trends in CDD, CWD, PRCPTOT, RIOmm, R20mm, R25mm 
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Figure 9. Time series of positive (negative) trends of four precipitation extremes in different stations with 
different climates. 



In many regions of the country, including regions 
in the northwest and west of the country, we were 
found negative trends in the consecutive wet days 
(CWD) index. The strongest significant negative 
trends were found in the stations of Gorgan (Caspian 
Sea) and Bam (arid and semiarid regions). In other 
regions, time series of the index show very weak 
positive trends were statistically non- significant. 
Almost every station in the west, north and middle of 
the country, especially in the west of Zagros 
Mountain and both side of Alborz Mountain, 
experienced very weak positive trends in the 



Consecutive Dry Days (CDD) index. In other words, 
there are not any significant trend and regular 
dispersal of CDD in these areas. The analysis 
indicated that the more stations in the south and 
southwest of the country experienced days with heavy 
precipitation in each year and positive trends for 
heavy precipitation were observed in these regions. 
There are a negative trend of heavy precipitation days 
(RIOmm) in the Northwestern, north and southeastern 
regions of the country. Furthermore, the highest 
number of heavy precipitation (RIOmm) was totally 
found at Tehran station in 1982. 
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Figure 10. Time series of positive (negative) trends of four precipitation extremes in different stations with 

different climates. 



The mega cities of the study area such as Tehran 
are different in many ways from other cities. The 
largest positive trends in R20mm have been found in 
Tehran station (2.44 days decade '). Possible reasons 
for having such patterns might be due to urbanization 
and/or increase of aerosols. In addition factors such as 
the local and regional air pollution and different 
urbanization characteristics may influence the nature 
and magnitude of the climatic trends in different 
ways. It is noteworthy that we observed significant 
positive trends in RIOmm, R20mm and R25mm, in 



the Bushehr station. This result suggests a more 
activity of rainfall systems such as Sudani cyclonic 
systems over these regions in recent years (Fig.8). 

It is important that the arid and semi-arid regions 
have experienced positive trends for maximum 1-day 
precipitation (RxlDAY) and this pattern might be due 
to drought several times especially in recent years. 

Positive trends for the RxlDAY index were 
observed in a vast area located in the north and west 
of the country (Fig. 10). The highest positive and 
significant trend was recorded in Esfahan (1.88 mm 
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per decade) located in the semi-arid region. The 
positive trends for maximum 5 -day precipitation 
(Rx5DAY) were found in Caspian Sea region and 
southern parts of the Alborz Mountain and southwest 
of the country. We found the steepest positive trends 
in Babolsar and Bushehr with slopes of 2.3 and 1.7 
mm per decade, respectively. On the other hand, 
negative trends for this index were found in the 
northwest and southeast regions of the country. 

3.3. Regional Extremes 

The regional averaged indices were computed as 
the mean of the indices at individual stations relative 
to their climatology. The significant trends in 18 
indices found in central part of the country (E region) 
are in accordance with the results of Evans et al. 
(2010) and the highest Variability of extreme 
temperature and precipitation was also found in this 
region. In the next stage, the largest variability 
occurred in Alborz and Zagros Mountain ranges and 
Azarbaijan respectively. Marked positive trends for 
indices like Warm nights (TN90p), tropical nights 
(TR20) and Warm Spell Duration Index (WSDI) were 
found over most regions of Iran. The largest negative 
trend of Warm nights (TN90p) of 1 percent per 
decade occurred in most regions of the country. In 
some regions like Azarbaijan and central part of the 
country time series of this index show the strongest 
positive trends. Generally, in parallel with Zhang et al. 
(2005) findings, on average, the Tropical nights 
(TR20) increased by 1 day per decade for the whole 
country. The strongest positive trends for this index 
have been found in central part of country and Zagros 
with trends of 1.82 and 0.85 day per decade, 
respectively. The whole country, on average, 
experienced an increasing trend of 1 day per decade 
for the Warm spell duration indicator (WSDI). 

We also found positive trends of GSL for the 
majority of the country, with the strongest in Zagros 
and Alborz Mountains and arid regions, showing an 
increase of 1.5 days per decade. The whole country, 
except some small areas in the Azarbaijan region, 
experienced a positive trend for Simple daily intensity 
index (SDII) and showed a trend of 2 mm per decade 
and is in accordance with the results of ShiftehSomee 
et al., (2012).We found positive trends of 1.1 days per 
decade for TMIN mean over the whole country 
except in Azarbaijan. In addition, the significant 
decrease of cool nights (TNI OP) was obtained over 
central parts and Alborz and southern regions 
respectively, at the rate of 1 .2 percent per decade. On 
the contrary, the magnitude of the positive trend in 
warm days (TX90P) found in this study for Zagros 



Mountain range (1.7 percent per decade), central parts 
and southern regions. The positive trends for TNx 
have been found in Zagros and Alborz mountain 
ranges and central parts with trend of 1.9 °C per 
decade. 

There were some areas in Zagros, Alborz and 
central part of the country that showed negative 
trends for diurnal temperature range (DTR). 
Maximum decrease in DTR was found in Zagros at 
the rate of 2.07 °C per decade. This decrease in DTR 
can be attributed to asymmetric rise in daily 
temperature extremes, i.e., TN is rising faster than the 
TX. Decreasing DTR found in Iran is consistent with 
the trend observed in different countries and at global 
scale (Karl et al., 1991. Decreasing DTR is related to 
fall in crop yield and therefore have severe 
implications for food production and food security 
both at regional and global scale (Peng et al., 2004; 
Peng et al., 2013; Rehmani et al., 2014). 

We observed a negative trend in CSDI over the 
Alborz mountain ranges and central parts, and 
negative trend in CWD over the Azarbaijan and 
Western regions of the country, although a clear 
negative trend was observed for ID0 and TX10P over 
central parts of the country. A positive trend was 
observed for the Alborz mountain ranges and central 
parts of the country for summer days (SU25) and a 
negative trend was observed for PRCPTOT over 
Azarbaijan regions. We found positive trends in 
TXmean for the central parts of the country. The 
results showed that most of the trends in the other 
indices were not significant over the whole country. 

4. Conclusions 

This study investigated long-term trends in the 
time series of the temperature and precipitation 
extreme indices at 30 synoptic stations in Iran on the 
daily timescales for the period 1961-2006. In this 
paper, the Mann-Kendall and the Durbin Watson tests 
were used to investigate the spatiotemporal trends. 
The findings of this study confirmed the results of 
previous investigation and present the new results for 
extreme temperature and precipitation in Iran. The 
positive trends in the temperature series is more 
consistent with the global warming events and have a 
significant effect on trends of extreme temperature 
events in recent years. The distribution of the trends 
for minimum temperature trends indicated that the 
negative trends mainly occurred in the most parts of 
Iran except in northwest part. On the contrary, 
positive trends were detected in more stations and in 
some areas, for example, from the central regions of 
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the country. No substantial differences were found in 
case of annual total precipitation and annual total 
precipitation from wet days (PRCPTOT), however 
substantial anomaly of such wet days occurred in 
1970. There are a large number and significant trends 
of heavy precipitation days in the south and 
southwest of the Iran. This result suggests a more 
increase of rainfall storms over south, southwest and 
central regions of the country. However, more 
statistically significant cases were found only at 
central part of the country and the highest number of 
stations with significant trends was also found in this 
region. Almost every station in the Azarbaijan region 
experienced less variability at regional scales. Further 
studies aimed to compare the trends found in our 
study at different spatiotemporal scales will be of 
great scientific interest. The results also suggest the 
need for further investigation on local anthropogenic 
intervention in the environment, which could be one 
of the major causes of climate change. 
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